Tunneling Effect Near Weakly Isolated Horizon 
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The tunneling effect near a weakly isolated horizon (Wflf) has been studied. By applying the null 
geodesic method of Parikh and Wilczek and Hamilton-Jacibi method of Angheben et al. to a weakly 
isolated horizon, we recover the semiclassical emission rate in the tunneling process. We show that 
the tunneling effect exists in a wide class of spacetimes admitting weakly isolated horizons. The 
general thermodynamic nature of WfH is then confirmed. 
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I. INTRODUCTION 



Since the discovery of the first exact solution of Ein- 
stein equations, studying properties of black hole is al- 
ways a highlight of gravitational physics. Laws of black 
hole mechanics and the famous Hawking radiation re- 
veal deep connections between classical general relativ- 
ity, quantum physics, and statistical mechanics. It is 
surprising that all stationary black hole can be described 
by several elegant laws. Furthermore, black holes are 
playing a major role in relativistic astrophysics, provid- 
ing mechanisms to fuel the most powerful engines in the 
cosmos. In experimental areas, black holes are also im- 
portant sources for gravitational wave detection. Unfor- 
tunately, the classical definition of black hole can not 
satisfy the requirements of practical research[l5. We 
need a "quasi-local" definition for black hole. In the 
last decade, much work has been done in this areafl*!- 
[l^. In contrast to event horizons, the quasi- local treat- 
ment doesn't require the global description of spacetime. 
Therefore, black hole mechanical laws can be general- 
ized to horizons while ADM quantities are replaced by 
local defined quantities. A natural question is whether 
these mechanical laws have thermal dynamical explana- 
tion. For stationary black holes, the answer is positive 
due to the discovery of Hawking radiation UM- Gen- 
eralization to isolated horizons is not straightforward be- 
cause Hawking's method requires the knowledge of the 
global geometry in a space-time, not just the geometry 
near a horizon. 

Recently, a new technique concerning Hawking 
radiation has been developed by Wilczek and his 
colleagues [l^ [13, [H, [l^. By treating Hawing radiation 
as a tunneling process near the Schwarzschild horizon and 
using WKB approximation, a semi-classical emission rate 



was derived as 
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where ASb~h is the change of the Bekenstein-Hawking 
entropy. Although this result is not exactly the thermal 
spectrum as first discovered by Hawking, it is consistent 
with the conservation of energy during the process. This 
method has been successfully extended to some other 
types of black holes [20l [2ll [2^ . All these calculations 
rely on specific forms of space-time. Since the tunnel- 
ing effect only concerns the physics around a horizon, it 
is natural to guess that tunneling effect could be a gen- 
eral phenomenon for all isolated horizons. The isolated 
horizon theory developed by Ashtekar et al. provides 
a perfect tool to explore this issue. By investigating the 
outgoing null geodesies near the horizon and applying the 
first law of the WIH, we recover Eq. ^ for the WIH. The 
WIH contains the essential features of a black hole, in- 
cluding mass, change and angular momentum. Thus, our 
work generalizes the result of the tunneling effect from 
stationary black holes to a wide class of spacetimes that 
admit WIHs. 

Angheben et al. reconsidered this problem using 
Hamilton- Jacobi ansatz. The Hamilton-Jacobi provides 
an alternative derivation of the tunneling effect. It is 
also a local method involving only an infinitesimal region 
around the horizon. Similarly, we apply this approach to 
a WIH, without the knowledge of the spacetime metric, 
we find the emission rate agreews with Eq. ([T]). There- 
fore, it confirms that the Hawking-like radiation exists in 
very general background. 

This paper is organized as follows. In section |TT1 
we first give a brief review of Ashtekar's weakly isolated 
horizon theory. Technical details can be found in ref . [l3| 
and references therein. In relation to WIH, we intro- 
duce the Bondi-like coordinates and calculate the outgo- 
ing null geodesic using Newman-Panrose formulism. Sec- 
tion [HTj contains the calculations on the tunneling effect. 
In section lTlI Al and llll Bl Wilczek's null geodesic method 
and Angheben's Hamilton-Jacobi method are employed 
respectively to study the tunneling effect near a WIH. 
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By using the result in section [TTl we find that both ap- 
proaches lead to the desired spectrum Eq. ((J). In sec- 
tion [11101 we discuss the role of coordinates in the cal- 
culations. The master equation in the Parikh-Wilczek 
method apparently depends on the choice of coor- 
dinates. The results in IIII Al apparently depend on the 
Bondi-like coordinate system chosen in section|lTl In sec- 
tion llll Cl we analyze the freedom of the coordinates that 
keeps the results invariant. We show that what is rele- 
vant in the calculation is not the coordinates, but the 
vector field approaching the horizon. Finally, we make 
some concluding remarks in section ITVl 



II. 



GEOMETRY OF WEAKLY ISOLATED 
HORIZON 



Base on works by Ashtekar and other authors [vj-p^. 
the weakly isolated horizon is defined by 

Definition 1 Let (M, g) he a space-time. Ti is a 3-dim 
null hyper-surface in M and Z° is the tangent vector filed 
of the generator ofH. Ti is said to be a weakly isolated 

h. oY\-Lon(WIH), if 

1) Tl has the topology of S"^ x 3?, 

2) The expansion of the null generator ofJi is zero, 

i. e. Qi — on H, 

3) TabV^ is future causal for any future causal vector 
and Einstein equation holds in a neighborhood ofH, 

4) [Cl , Da]l'' — on TC, where Da is the induced 
covariant derivative on Ti.. 

Since we are going to study the null geodesic behavior 
around the horizon, we need non-singular coordinates 
to describe the geometry through the horizon. For the 
Schwarzschild black hole, Parikh and Wilczek chose 
the Painleve coordinates. For the WIH, a convenient 
choice is the Bondi-like coordinates [2^. Such coordi- 
nates have been used to prove the local existence of WIH 
[2^ . The coordinates are constructed as following: The 
tangent vector of null generator of H is 1°". Another real 
null vector field is n° . The foliation of Ti. gives us the nat- 
ural coordinates {9, (f). Using a parameter u of P and Lie 
dragging (6*, 0) along each generator of Ti, we have coor- 
dinates (m, 0, 0) on Ti. Choose the affine parameter r of 
n" as the forth coordinate, we get the Bondi-like coordi- 
nates {u, r, 9, (j)) near the horizon. In these coordinates, 
we choose the null tetrad as 



r = 



du dr dQ d( 

a 



m 



dr 



d 



d 



dr dC, dQ ' 



(2) 



where C/ = X= w = 0on7i (following the notation 
in ref.[l3|, equalities restricted to Ti will be denoted by 
" = " ) I C = e*"^ cot I . In words of Newman-Penrose for- 
mulism, we also require the above tetrad satisfy following 
gauge, 

v = ^j = T = a + j3 — TT^ji — fi, — 

e — £ = K = 0, (3) 

which means these tetrad vectors are parallelly trans- 
ported along n° in space-time and Lie drag-free along 1°' 
on Ti. Furthermore, the forth requirement in the def- 
inition of WIH implies there exists a one form oja on 
Ti such that Dal^= ujal^ and Ciuoa^ 0. In terms of 
the Newman-Penrose formulism, LOa can be expressed as 
uja = + £)na + (a + /3)ma + (a + I3)ma. The above 
equation means {£ + £) is constant on Ti. In such a space- 
time, the out-going null geodesic is described by We 
have known that U =0. Commutative relation of /° and 
n° tells us that 



dr 



(e + e) 



TTLO, 



(4) 



which means Hr-^(s - 

or ^ 

U near Ti is 



e). Then the behavior of function 



U — {e e)r -f o(r). 



(5) 



Based on the discussion in ref.[13|, not any choice of 
time direction can give a Hamiltonian evolution. In other 
words, only some suitably chosen time direction can lead 
to a well-defined horizon mass. In ref.[3], Ashtekar et 
al. gave a canonical way to choose the time direction t"' 
for a WIH. Compared with the Schwarzschild case 
the parameter of takes the place of the Killing time. 
Using the tetrad in Eq. ([2]), We find the time derivative 
of r along the outgoing geodesic 



du dr 
dt du 



o{l))U^Bt{e- 



- e)r - 



o{r). (6) 



By definition, the surface gravity oiTi is nt ■— Btl°'uJa = 
Bt{e-\-e). Because Bt{e-\-s) is constant on Ti, the zeroth 
law of black hole mechanics is valid for WIH. 

With the canonical time direction f^, Ashtekar et 
a/.jl3j showed that the following relation holds on Ti. 



SM. 



it) _ lit 



n 



(7) 



dc dc 



where M^^-* is the horizon mass, a-^ is the area of the cross 
section of WIH, fit is the angular velocity of the horizon 
and J-H = —-^§g{uJa'il^°')dS is the angular momentum. 
This is the first law of black hole mechanics for WIH, 
which is a generalization of the first law for stationary 
black holes. 
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III. TUNNELING EFFECT NEAR WEAKLY 
ISOLATED HORIZON 



A. Null geodesic method 

In this section, we consider the tunnehng effect 
near a weakly isolated horizon. We apply Parikh and 
Wilczek's semi-classical method [H [Isl, [l^ to a WIH. 
We shall use the Bondi-like coordinates introduced in 
section HIl It has been demonstrated [1 7] that the WKB 
approximation is justified near a horizon. From the WKB 
formula, the emission rate F can be expressed as 



F ~ exp(-2/mS'). 



(8) 



If the positive energy particle goes outwards, from ri„ to 
fouti the imaginary part of action is given by 



ImS — Im 
= Im 



Prdr, 



dprdr. 



(9) 



From the Hamiltonian equation, pr can be expressed as 

dH 



dpr ~ 



With the help of Eq. dTO]) and Eq. 
of the action becomes 



(10) 

the imaginary part 



ImS 



Im 



Im 



M 

M- 



dH 



dr 



M 



dr 



Bt{e 



- e)r 



o{r) 



dH 



dH 



Bt{e + e) 



(11) 

Following the argument in [16|, , we fix the total mass 
of the space-time and vary the black hole mass. In section 
im the general form of the first law is given by Eq.©. For 
the spherically symmetric case (J7t = 0), we can express 
dH in terms of the surface gravity and the horizon 
area a-n, 



dH = -du!^^ = ~— dan 



The minus sign in above equation comes from the con- 
servation of energy. Then we have 



ImS 



dH 



Btie + e) 



d (—) 

2TrBt{e + e) V 4 / 



(13) 



and the emission rate F is 



F - exp(-2/mS') = exp(A^) = exp(AS'i3_H) . (14) 

This gives the thermal spectrum of radiation, which is in 
agreement with Eq. ([T|). 

For an axial symmetric horizon, the first law 
becomes [13, [H 



dM-H = ^ dan + dJn- 



(15) 



Using the formula in ref . [22] and the first law ([7]) , we have 



ImS = Im j [L- P^^)dt 




(16) 



where the conservation of angular momentum has been 
used in the fifth equality. Eq. (flB)) gives the same emis- 
sion rate near a rotating horizon. Therefore, we have 
proven that the semiclassical emission rate in form of 
Eq. ([T]) holds for a general WIH. 



B. Hamilton-Jacobi method 



In this section, we use the Hamilton-Jacobi method 
to re-calculate the Hawking radiation near a WIH. This 
method has been used in stationary space-times pol. [2]|. 
We need first modify Bondi-like coordinates (i, r, 6*, (/>). In 
general, we can always write t° as t°'=Btl°' — Xq ( gfr) , 
where_A5 = 3,4, = BtX\n, X^ = BtX\n, x^'= (, 
— (. The modified Bondi coordinates are constructed 
in following way: Instead of using u, we use parameter 
of t"' as the time coordinate; the coordinates {6, (p) are 
still obtained from a foliation of Ti and extended to H 
by Lie-dragged along t° ; coordinate r is chosen such that 
n° = / (^) . Under these coordinates, the only different 
gauge is that X is constant on H. In these coordinates, 
we write the tetrad as 



J_d_ 
B'tdi 
d 



4-4 

or 



X 



d_ 

5C' 



/ 



dr ' 



4 



d ,o 9 ^. d 



d 



d 



d 



(17) 



where BtX=Xo, f=l, uj^Q, U=0. Then the metric g^"' 
takes the form 







Bt 







Bt 



fx^-{Lj^^+uj^^) -(c-^e^ + e^c^^ 



(18) 



In addition, we shall use the same gauge condition ([3]). 

Based on this ansatz, the action function / of an 
outgoing particle should satisfy following equation : 

,2 



0. 



(19) 



Comparing with null geodesic method, we also need 
to control other components of metric near the horizon. 
Note that the gauge choice in section HI] has fixed the 
metric components on Ti.. The behavior of metric compo- 
nents near the horizon is controlled by Cartan structure 
equations. 



dU_ 

dr 
dX 
dr 
du! 
dr 

de 

dr 
dr 



+ e)- 


TTLO 






_ f 






flCU 




'T 






f 








/ ' 



Dlnf, 



(20) 



where D :— l°'Va, 5 := m°-V a- Here we have used the 
modified Bondi coordinates {t,r,9,(l)). The behavior of 
metric near horizon is 



U — {e + e)r + o{r), 
X^Xo + 0(r), 

e' = o(i), 

iv = 0{r). 

The time derivatives of these quantities are 

dU 

'dt 

dX 

'dt ' 
duj 
'dt ' 
dX"^ 



(21) 



dt 



= o(r), 
0(r), 
0(r), 
--0{r). 



(22) 



We consider the variable separation solution for / 



I = V{t) + W{r) + Jix"^). 
Eg. pQ)) then becomes 

2{fU-\uo\''){W'f 

+2[(/X^ - oj^ - Loi^)JA + -^VW 

-{e~e+~ee)jAJB+m^ 



(23) 



Bt 
0, 



(24) 



where " • " means 
Consequently, we have 

W{r) = 

where 



means 7^- and J a 



dJ 



-B + ^/B'^ ~ AC 
'A 



dr, 



(25) 



A = 2(/C/-|c.n, 



/ 



B = ifX^ - Cue - ^nJA + jrV^ 

Bf 

C = _(e^P+e^^'')JAJB+m2. 
Taking the time derivative on both sides of Eq. (P5)) gives 

2{fu + fir -QCo- uju)){w'f 

+2 \{jX^ + fX^ - ~ - iu^ - u;^^)Ja 



f ■ f ■■ 

'—V + —V 
Bt Bt 



w 



(26) 



Because Eq. ([26)) holds smoothly in the neighborhood of 
H., with the help of the asymptotic behavior of the tetrad 
(1^ and (HH), we get X^Ja + I; = 0. Furthermore, 
similar derivation shows dsiX^JA) = 0, i.e. X'^Ja + 
-g^=const. In the axisymmetric case, Xq is constant on 
7i[13]. Thus, V and J a are constants. Comparing with 
the stationary case, E = — is the energy of this particle 
and Ja is the angular momentum of the particle. The 
energy condition guarantees E — BtX^JA > 0. With the 
above results and the asymptotic behavior, we find that 
integral (j25p has a simple pole at horizon. The imaginary 
part of the action function is determined by this pole. 
Specifically, 



Im{I) ^ Iml W'dr = n 



E — BiXqJa 



(27) 



where K( = Bt{e + e). Combining with the WKB as- 
sumption ^'(x) — exp(j/), we get the thermal spectrum 
r ^ exp[— /3(i? — BiX^Ja)]- As we emphasized in sec- 
tion ini Only a canonical time direction t° can lead to 
a well-defined horizon mass and angular momentum, as 
well as the black hole mechanical law. For such an ob- 
server, t°'=Btl°' — r2t'0°! where -0" is a Killing vector of 
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Ti and tangent to the leaf S, fit is the angular velocity of 
the horizon. In this case, we choose Xq as = flt^^- 
Because of the conservation of total energy and angular 
momentum, E is the variation of the energy of black hole 
and ?A(J) is the variation of the black hole. Combining 
with the quasi- local black hole mechanical law ([7]), we 
obtain T - exp(A2|t). 



Note that gtt=0 and the metric is non-degenerate around 
the horizon. Thus gtr must be non-zero at the horizon. 
So we can simplify Eq. (|29p by Taylor expanding the 
square root term around g^^.. Then wc find 



2.9tr 



gABX'^i^ 



gtt 



(31) 



C. Freedom of coordinates 

The calculation in section IIII Al appears to rely on 
a specific coordinate system. The purpose of this sub- 
section is to investigate how much the result depends on 
the choice of coordinates. In order to calculate the tun- 
nelling effect, the crucial step is to calculate the residue 
of 1/r. We consider the coordinates {t,r,9,(f>) as defined 
in section [ill where t is the canonical time. It will be 
sufficient to focus on the following ine element 

ds^ = gttdt^ + 2gtrdtdr + grrdr^ + gAsdx^dx^ . (28) 

These coordinates are nonsingular through the horizon. 
By letting ds^ — 0, f, the derivative of r along an outgo- 
ing null geodesic, can be easily solved as 

• _ -gtr + Vgtr ~ 9rr{gtt + gABX^X^ , . 



Since (^)" approaches the horizon, gtt must vanish on 
the horizon. Since r is chosen to be a constant on the 
horizon, r also vanishes on the horizon. Then we have 

grrigtt + gABx'^x^)=0 ■ (30) 



It is not difficult to see that the approximation made in 
Eq. ([3T|l preserves the residue. Then Eq. (f3T|) leads to 
the desired emission rate as shown in section IIII Al Now 
consider the coordinate transformation 

t = T + fiR) 

r = R, (32) 

where /(i?) is a smooth function around the horizon. 
Under the new coordinates, the metric becomes 

ds^ = gttdT^ + [2gtr + gttfir)]dTdR 

+ [grr + 2gtrf'{r)+gttnrf]dR^ 
+gABdx'^dx^, (33) 
which gives the radial null geodesic equation 

[grr + 2gtrf'{r)+gttf'{rf]R^ 

+ [2gtr + gttf{r)]R + gtt + gABX^x^ = , (34) 

where R = dR/dT. The outgoing solution for -i is 



1 
R 



2{grr + 2gtrf{r)+gtt.nry 



-2gtr - gttf'ir) + y/i2gtr + gttf'irW ~ ^9rr + 2gtrf'{r) + ff«/'(r)2] [gtt + gABX^x^] 



(35) 



Note that gtt = at the horizon. Since the coordinates 
are nonsingular, gtr cannot be zero at the horizon. In 
order to find the residue of l/R, we Taylor expand the 
square root around (2gtr + gttf'ir))'^ to order one (terms 
with higher orders are not relevant to the calculation of 
residue). Then Eq. ((35)) becomes 



1 

R 



2gtr + guf'{r) 
gtt + gABX^i^ 



(36) 



Noting that 1744 vanishes on the horizon and comparing 
Eq. (|36p with Eq. (|3ip. we see that 1/R has the same 
residue as 1/r. 



Now we are going to show that the result is invariant 
under a rescaling of r. Consider the following coordinate 
transformation: 



t 



h{f) , 



(37) 



where f is any smooth function of r satisfying f = on 
the horizon. Similarly, we find 



f7 + Vte 
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9tt + QABX^i^ ' 

Suppose the residue of 1/r at the horizon is a. Then 
1 a/'(r) af{r) af'{f = 0) a 



(38) 



/(f) /'(f = 0)f f 



(39) 



which shows that 1/r and 1/f have the same residue on 
the horizon. 

We have shown that there exist two classes of co- 
ordinate transformations, Eq. and Eq. ([57)1 . that 
keep the residue invariant in the tunneUing effect. How- 
ever, the two transformations share a common feature: 
they both preserve the tangent field (^) . So we may 
conclude that the emission rate is determined only by 
a specific family of orbits approaching the horizon and 
parameterized by some coordinate time t (which is the 
Killing time in the Schwarzschild case) . All comoving co- 
ordinate systems preserving the tangent field (^)" are 
equivalent in the sense that they lead to the same emis- 
sion rate. 

Now we keep the orbits fixed on the horizon but 
vary the orbits off the horizon. This can be realized by 
the following coordinate transformation 



r = h{T,R) 
t = T. 



where A = |^=0. Then 



d_ 

df 



dr 



(40) 



(41) 



which shows that (^)° and (^)° coincide on the hori- 
zon but are different off the horizon. By a similar deriva- 
tion, we have 



1 

R gtr + 9rrh - Vsir ~ 9rr{gu + gABX'^X^) 

_ -2gtrh' 

2gtrh + gtt + gABX^i^ 



(42) 



A derivation similar to Eq. ([39|) shows that h' in Eq. 
(|^^ docs not change the residue. However, if h ^ r near 
the horizon, it will cause a change of residue. For h ^ 
[n > 2), the residue keeps invariant. Since h measures 
the difference of the two families of orbits approaching 
the horizon (see Eq. (|4T|l ). the above results show that 
when the difference is in the order of r, the residue will 
change. But if the difference is smaller than r (higher 
order of r), the residue is unchanged. Thus, if the two 
families of orbits in Eq. (|¥T|l satisfy h , n > 2, we 



call them equivalent in the sense that they lead to the 
same residue. 

IV. DISCUSSION 

The Hawking-like radiation has been derived as a 
tunneling process near weakly isolated horizons. The null 
geodesic method and Hamilton-Jacobi ansatz lead to the 
same result. This indicates that thermal radiation is a 
generic property of horizon, not only for stationary black 
holes. The choice of time direction has played a key role 
in both methods. However, there are subtle differences 
in the two methods. In the null geodesic method, only 
the canonical time direction can define the horizon mass 
and lead to the first law of black hole mechanics. In the 
Hamilton-Jacobi method, the thermal spectrum exists for 
any choice of time direction. However, only when the 
canonical time direction is chosen, can the Kt term in the 
expression be interpreted as the Hawking temperature. 
We also find that different observers will give different 
temperatures even if they have the same limit curve at 
the horizon. We find the difference dt — dt^ ~ O(r^) can 
insure the invariance of the temperature, where tc is the 
canonical time given by Ashtekar et al.^^. This result 
is reasonable because the energy and angular momen- 
tum measured by arbitrary observers cannot be directly 
related to the horizon mass. The canonical time direc- 
tion corresponds to the Killing observers in Schwarzschild 
spacetime, where other observers are like the Rindlcr ob- 
servers of the Unruh effect. 

Our discussion in this paper is confined to vacuum 
solutions. There is no obvious difhculty to extend the dis- 
cussion to matter fields, for example, the Maxwell field or 
other gauge fields. When matter is considered, the field 
equations can be controlled in similar ways and the first 
law including mater fields is also known [l|-[l3j. The next 
interesting thing is studying dynamical horizons. Some 
special cases have been investigated ,25.]. The general cal- 
culation on dynamical horizons will be done in future 
work. 
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